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ABSTRACT 

The properties of tidally induced arms provide a means to study molecular cloud formation and 
the subsequent star formation under environmental conditions which in principle are different from 
quasi stationary spiral arms. We report the properties of a newly discovered molecular gas arm 
of likely tidal origin at the south of NGC 4039 and the overlap region in the Antennae galaxies, 
with a resolution of 1"68 x 0'.'85, using the Atacama Large Millimeter/submillimeter Array science 
verification CO(2-l) data. The arm extends 3.4 kpc (34") and is characterized by widths of < 200 pc 
(2") and velocity widths of typically AV ~ 10-20 km s . About 10 clumps are strung out along this 
structure, most of them unresolved, with average surface densities of E gas ~ 10-100 M Q pc -2 , and 
masses of (1-8) xlO 6 Mq. These structures resemble the morphology of beads on a string, with an 
almost equidistant separation between the beads of about 350 pc, which may represent a characteristic 
separation scale for giant molecular associations. We find that the star formation efficiency at a 
resolution of 6" (600 pc) is in general a factor of 10 higher than in disk galaxies and other tidal 
arms and bridges. This arm is linked, based on the distribution and kinematics, to the base of the 
western spiral arm of NGC 4039, but its morphology is different to that predicted by high-resolution 
simulations of the Antennae galaxies. 

Subject headings: galaxies: individual (Antennae galaxies) - galaxies: ISM - galaxies: structure - 
stars: formation 



1. INTRODUCTION 

The gas and star formation (SF) properties in arms 
induced by tidal interactions are believed to differ from 
quas i-stationary s piral arms in spiral density wave theory 
(e.g. Dobbs|20lT ). However, the studies of the molecular 
properties in tidal arms are quite limited, even though 
they are important elements to shed light into our ideas 
of cloud fragmentation and SF process. In addition, the 
formation of these features contributes to the increase 
of star formation ra te (SFR) in interacting galaxies (e.g. 



Teyssier et al. 2010). Existing studies are focused on tidal 
tails, where the molecular g traced by CO, is only 
detected in a few rather disconnected clumps, and usu- 



ally within H I tails (e.g. | Tayloi et al.||2001 



al. 



2004 



Lisenfeld et 
re SI 1 ' prop- 

erties" estimated in several collisional debris show similar 



Walter et al. 200t>p. Surprisingly, t 



values to those measured in galactic disks (praine et al. 



2001 Boquien et al. 



20TT||Duc fc Renaud||201l| . 



Une ot the best studied colliding systems is the An- 
tennae galaxy (NGC 4038/9, Figure [l]) , partly because 



it is o nly at 22 Mpc distance from us (|Schweizer et al 
2008 1 , and thus it is a good target to study the properties 
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of molecular gas in tidal arms with high spatial resolu- 
tion (1" roughly corresponds to 100 pc). The two spiral 
gala xies started to interact only 400-600 million years 
ago ( Mihos et al.|1993 Teyssier et al.|2010 1 , making this 
system one of the youngest examples ot an on-going ma- 
jor galaxy merger. Two prominent H I tidal tails extend 



about 70 kpc north and south (e.g. Hibbard et aL~|200T ) 
which have been successfully used among oth er informa- 
tio n to c onstrain numerical simulations (e.g. Teyssier et 
aO20M. 

High-resolution maps (< 5") of the thr ee lowest tran- 
sitions of CO have been prese nted by Wilson et al. 
( |2000| , |Ueda et al.| ( |2012| > and |Wei et al.| Q2012p us- 
ing the Uwens Valley Radio Observatory (OVRU), the 
Plateau de Bure Interferometer (PdBI), and the Submil- 
limeter Array (SMA). The molecular gas as traced by 
CO emission is mostly distributed around the two nuclei 
of NGC 4038/9, a region to the west of NGC 4038, and 
along an overlap region between them, and co exists with 
a large amount of super stellar clusters (SSCs Whitmore 
& Schweizer||1995|). 



In this Letter, we report the molecular gas properties 
of a tidally induced molecular arm found in NGC 4039 
that had not been previously detected. We have selected 
this region because it is a quite narrow and elongated 
molecular feature that allows us to study giant molecular 
cloud formation in a filamentary structure. In addition, 
it is not affected by projection with other components. 

2. CO(2-l) ALMA OBSERVATIONS AND DATA REDUCTION 

We use the Atacama Large Millimeter/submillimeter 
Array (ALMA) science verification CO(2-l) data of the 
Antennae galaxies, retrieved from the ALMA science 
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portaj^] Band 6 observations with 14 antennas (11 North 
American and 3 East Asian antennas) were carried out 
in 2011 May and June, using two spectral windows of 
1.875 GHz bandwidth (or -1600 km s" 1 ) and 3840 0.488 
MHz channels, and with one of the spectral windows cen- 
tered on the CO(2-l) line. The observations were split to 
cover separately the northern and southern galaxy, with 
13 and 17 pointings, respectively. In this Letter, we fo- 
cus on the southern mosaic, which includes the overlap 
region and NGC 4039. The synthesized beam is charac- 
terized by a FWHM size of l'/68 x 0"85 and a position 
angle of 78°. The covered area was approximately 1.'5 
x Of 9, along a position angle of 60° (Figure [I]). The 
observed UV range was 14 - 236 m. 

Editing, calibration, and imaging were performed in a 
standard manner using CASAPj similar to that presented 
in the CASA guide[^]for the CO(3— 2) science verification 
data of Antennae galaxie^] Due to the large complex- 
ity of the emission we restricted the cleaning to manu- 
ally selected boxes where emission was apparent for each 
channel. The rms of a free-line 20 km s^ 1 channel is 
1.2 mjy beam -1 . However, this image is dynamic range 
limited due to the bright emission in NGC 4039 and the 
overlap region, and we have estimated that the rms in 
each channel close to these features is — 6mJybeam _1 . 
The absolute flux accuracy is estimated to be about 10%. 

Figure ^1 shows the integrated intensity (moment 0) , 
velocity field (moment 1), and velocity dispersion (mo- 
ment 2) maps. To separate real emission from noise, 
we computed only those regions which show emission in 
three consecutive channels of 10 km s" 1 in a data cube 
that was convolved to 5 pixels (Of! 2 pixel -1 ). The veloc- 
ity range used is 1300 - 1800 km s" 1 . In addition, we 
used a clipping of 1.2 mjy beam" 1 per channel for the 
moment map, and 1.5 mjy beam -1 for the moment 1 
and 2 maps. 

3. DISTRIBUTION AND KINEMATICS OF THE 
MOLECULAR ARM 

Figure [2] shows that the molecular emission is partic- 
ularly bright in the NGC 4039 nucleus and the overlap 
region between the two galaxies, where it appears con- 
centrated in five large mo lecular complexes, in agreement 
with previous CO map s ( Wilson et al.|2000| Ueda et al. 
2012 Wei et al. 2012). A comparison with these pre- 
vious maps, as well as ALMA CO(3-2) data, showed 
that the brightest complexes are well collocated. The 
total CO(2-l) flux of the entire southern mosaic map is 
1460 Jy km s . The total flux of the map is s lightly 



larger than previous estimates by Wei et al. (2012) 
1358 Jy km s _1 . By comparing with single dish obser- 
vations, |Wei et aL"1 infer that flux loss in their interfer- 
ometric observations is —12%, which is comparable to 
absolute total flux uncertainties. 

The CO(2-l) map (Figure [2]) shows a molecular arm 
in the southern region close to NGC 4039 and the 
overlap region. The arm extends 34" end to end, or 
3.4 kpc. Some regions of this feature are slightly re- 
solved in the N-S direction, but the width is typically 



6 http:/ /almascience. nao.ac.jp/almadata/sciver/AntennaeBand6/ 

7 http://casa.nrao.edu 

8 http:/ /casaguides.nrao.edu/index.php?title=AntennaeBand7 

9 http:/ /almascience. nao.ac.jp/almadata/sciver/AntennaeBand7/ 



unresolved at —200 pc. The CO(2-l) flux of this feature 
is 60 Jy km s" 1 , a 4% of the total flux in the ALMA 
CO(2-l) southern mosaic. This component is also par- 
tially seen in the ALMA CO (3-2) maps, although with 
a considerably lower signal-to-noise ratio. This may be 
due partly to low excitation, especially in regions further 
from NGC 4039 and the overlap region (Figure [l]). 

The velocity of this component ranges from 1600 to 
1750 kms -1 , but the velocity dispersion is typically 
—10-20 km s" 1 (Figure [2] lower left panel). This is 
slightly higher than the typical velocity dispersion of 5— 
10 k m s -1 in the molecular gas o f face-on galaxy disks 
(e.g. Combes & Becquaert 1997). The overlap region 
shows much larger velocity gradients and velo city dis- 
persions, up to —70 km s , in agreement with Ueda et 
af] ( |2012[ ) and |Wei et al.| ( |2012h . Figure [3] (upper panel) 
shows that the molecular arm s velocity centroids vary 
continuously and slowly within the velocity range 1650 
- 1750 km s _1 . This continuity is kept along the spiral 
arm of NGC 4039, but not in the vicinity of the brightest 
giant molecular gas complex in the overlap region, where 
several spectral components are found in projection. 

The molecular gas surface densities were derived 
using Sh2 (Mq pc -2 ) = 156 
10 20 ) S CO (2-i)Ay [Jy km s -1 arcse 
S CO (2-i)AU term is the CO(2-l) integrated flux den- 
sity. We adopted an intensity CO(2-l) to CO(1-0) ratio 
7(2 1)/I(1 0) = 1.0 and a CO-to- H 2 conversion fac - 



x (X/3.0 x 
2 1, where the 



tor X = 3.0 x 10 20 cm -2 (K km s" 1 )" 1 (Wei et al. 2012) 



We did not correct by 1.36 to account for He. Along the 
molecular arm, the gas surface densities of the brightest 
components are typically — 10-100 Mq pc -2 . Note that 
in the overlap region they are as high as 1000 Mq pc~ 2 . 

At least 10 clumps are seen along this structure as 
shown in the lower panel of Figure [3j The masses of these 
clumps span from (1-8) x 10 6 Mq, with the most mas- 
sive ones being closer to NGC 403 9. Using also CO (2- 
1) observations, Wei et al. (2012) identified a distinct 
break in the mass function of molecular clouds at — 
3 x 10 6 Mq, with higher values corresponding to ac- 
tively forming regions. The clumps in the molecular arm 
are therefore in between these two mass populations. For 
a cloud size scale of —100 pc, using a typical velocity dis- 
persion of 10 km s -1 , and a geometric factor a g for the 
density profile of 5/3 as in Wei et al. (2012), the virial 
mass is 7xl0 6 Mq. This suggests that while the molecu- 
lar clumps close to NGC 4039 are compatible with struc- 
tures that are bound gravitationally, we cannot discard 
that they are not fully bound further away. Higher res- 
olution observations to resolve these complexes are nec- 
essary to set stronger constraints. 

The separation between these complexes (peak to peak 
average value) is found to be almost equidistant at about 
350 pc, assuming that the effect of inclination is negligi- 
ble. It is remarkable the resemblance to classical "beads 



on a string" (e.g., Elmegreen & Elmegreen 1983) along 
the molecular arm. The characteristic length (A) to fila- 
ment diameter (D) ratio, X/D, is X/D — 2, assuming A 
= 350 pc and D = 200 pc. 

Although usually applied to different scales than those 
presented here, a number of models have been created to 
describe fragmentation of filamentary structures. If frag- 
mentation is caused by self-gravity the distance between 
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adjacent clumps is of the order of the Jeans length. The 
typical Jeans length depends on gas pressure (both ther- 
mal and turbulence), magnetic field and rotation. The 
expected X/D ratio obtained in these models, X/D, is 
typically 2-4, with smaller ratios found f or faster rota- 



tion and/ or stronger magnetic fields (e.g., Nakamura et 
aLl |1993| |Hanawa et al.j|1993| ), in agreement with the 



UO(2-lj observations. 

The filamentary structure of the molecular arm is 
unique to probe the characteristic separation between 
giant molecular associations (GMAs), and allows us to 
study the typical size scales for fragmentation occurring 
between parsec and kiloparsec scales. Observations of 
Galactic filamentary clouds show characteristic separa- 
tions of ab out three times their diameter (a few parsec 
scale; e.g., Schneider & Elmegreen 19791. The separa- 



tions of GMAs in nearby galaxies and our Galaxy are 
characterized by a large scatter, with separations of a 
few 100 pc to 1 kpc (e.g., [Rh^[l993l|Koda et al.|[2006|. 
These are located within larger scale complexes of SF 
regions separated on average a few kiloparsecs and typi- 
cal diameters of 0.5-1.0 kpc, in a configuratio n that also 
resembles beads on string morpholo gy (e.g. |Elmegreen| 
fc Elmegreen||1983| |Elmegreen et aT] 2006| |Smith et al. 
2010| |Mullan et al.|2 01ip. 



The general properties of the molecular arm are com- 
parable to that of spiral arms in non-interacting galaxies. 
However, the width seems smaller than the typical width 
in spiral arms, usually in the range of 400- 1500 parsec 
based on measu rements using Ha maps (e.g., [Kennicutt| 



& Hodge 1982). Although the inclination of this struc- 
ture is rather uncertain, we find it unlikely that the ob- 
served molecular arm can be reproduced with a single 
pitch angle, which may indicate that it has been created 
or affected by tidal interactions. 

4. STAR FORMATION LAW AND STELLAR CLUSTERS 

In order to investigate the SF properties in the molec- 
ular arm, we derived the SFR surface densities (£ SF r) 
from available FUV and 24 /im maps. We used FUV 



da ta fro m the Nearby Galaxy Survey (Gil de Paz et 
al. 2007). The FUV band covers the wavelength range 
1350-1750 A, and provides a resolution of 470. We cor- 
rected the FUV maps for f oreground extinctio n using 
A FUV = 8.24x E(B - -V) (|Wyder et all|2007], where 



E(B - — V) = 0.046 mag as estimated by|Schlegel et al. 



11998). We use the 24 /im map presented in Zhang et 
al.|(| 2010| ), which was obtained with the MIPS ( fHleke et| 
004 ) on board the Spitzer Space Telescope, and was 



corrected for background, matching of individual data 
frames, cosmic-ray removal, flat-fielding and mosaicking 
using MOPEX. The FWHM of the MIPS maps point- 
spread function (PSF) at 24 [im is 6". The FUV map was 
convolved to the same sp atial resolution of 6 ". The Esfr 
was calculated following Leroy et al. (2008): Esfr [M© 
yr -1 kpc -2 ] = 3.2 x 10 * hi [MJy str^] + 8.1 x 10 -2 
Ifuv [MJy str -1 ]. The sensitivity limit of the Esfr is 1 
x 10 -4 M© yr -1 kpc" 2 . 

We also convolved the molecular gas surface densities 
to 6" and converted the molecular gas and SFR surface 
density maps to a pixel size of 3". For the SFR map, a 3a 
threshold was used, or 0.018 M Q yr -1 . Figure [4 shows 
the pixel-to-pixel SF law plot for different regions in the 



southern mosaic. The average star formation efficiency 
(SFE), the SFR per unit of gas, for the molecular arm 
is log(SFE [yr -1 ]) = -8.2 ± 0.3, similar to that of the 
NGC 4039 region log(SFE [yr -1 ]) = -8.8 ± 0.3, and both 
are a factor of 10 lower than in the overlap region. These 
are still larger than the av erage SFE for normal galaxies 
of log(SFE [yr" 1 ]) = -9.3 ( |Komugi et aL|2005| [Bigiel et 
a O2008"l ). 

Note that for the molecular arm we used only regions 
in the SFR map that are more than 11" away from the 
brightest points in the overlap region and NGC 4039, in 
order to avoid confusion with the airy ring at ^8" in the 
24 fj,m PSF. The molecular arm region pixel closest to 
the overlap peak (611 MJy str -1 , with background re- 
moved) is about 13" away from the overlap peak. The 
24 jum PSF is about 1% here, or about 6 MJy str -1 . The 
actual 24 /im flux density at this region is ^24 MJy str -1 , 
so the PSF contamination in this case is about 25%, or 
about 0.1 dex on the SF law plot. Considering the rms 
dispersion of the SFE (0.27-0.46 dex), this can be con- 
sidered negligible. 

A large number of bright young star cl usters (^700 
identified) exist in the Antennae galax ies ( |Whitmore fc 

R 



Schweizer 1995 Whitmore et al 



1999) 



lowever, it is 

remarkable that the SSCs number density is quite low 
along the molecular arm. We can only identify five can- 
didates of SSCs within 5" of the molecular arm (see Fig- 
ure |3b - Based on V — I colors and th e color-age diagram 

(|2010 " [), we could only con - 



presented in Whitmore et al. 



litmore fc Schweizer|1995 ) 



strain the age for SSC 33 (W 

<10 Myr. Based on the relative location and young age 
of SSC 33, we suggest that it is likely associated with 
some of the more massive molecular gas complexes in 
the molecular arm. 

5. DISCUSSION: ORIGIN OF THE MOLECULAR ARM 

That the molecular arm is a single entity linked to 
the spiral arm of NGC 4039 is proven by the smooth 
and continuous velocity change along its different com- 
ponents in the diagram shown in Figure [3] Due to the 
link in terms of distribution and kinematics with the spi- 
ral arm of NGC 4039, we believe that the most plausible 
explanation is that the molecular arm is of tidal origin, 
presumably caused by the recent passage of the north- 
ern galaxy, NGC 4038. |\\ liit more et al.| ( |2010[ ) identified 
in Hubble Space Telescope (HST) images the molecular 
arm as a dark, relatively circular dust cloud (found near 
the so-called regions " Outer 3" and " Outer 6" in Figures 
20(c) and (g), and suggested that this material may be 
falling back into the galaxy. 

If the stripped material contained both atomic gas 
(H I) as well as molecular gas, then we would expect 
that H I is also extended in this region, in disagreement 
with the observations. The H I integrated density dis- 
tribution (Hibbard et al. 112001 1 does not show any clear 



correlation with the molecular arm, although the avail- 
able H I map's resolution is a factor of 10 worse than that 
of the CO(2-l) map. In the region of the molecular arm, 
the H I flux density is Shi < 0.03 Jy beam -1 kms -1 , or 
H I surface de nsity of lM^pc" 2 , following for exam- 



ple equation in Espada et al. (2011). This suggests that 
mostly molecular gas material was dragged. If not of 
tidal origin, another explanation might be a bubble ere- 
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ated by massive SF. However, this is unlikely due to the 
necessary ene rgy output, E ~ 10 51 J following a similar 
method as in Sakamoto ct al. 2006 > which would require 
~ 10 7 supernovae explosions. 

Next we compare with available high resolution nu- 
merical simulations of the Antennae galaxies. |Karl et ah] 
( 2010 ) recently presented simulations with resolutions of 
~70 pc, using radiative cooling, SF, feedback from su- 
pernovae, and a refined merger orbit. These simulations 
reproduce in general the observed gas morphology and 
kinematics, as well as the SFR. However, note that the 
predicted projected gas surface density map (their Fig- 
ure 5.3) shows significantl y more gas in the southern disk 
than in the H I maps, and Karl et aL] ascribe it to most of 
the gas being in molecular rather than in atomic phase. 
The ALMA maps show such missing molecular gas, but 
the morphology is different. We notice that there is a fil- 
amentary structure close to the east of the overlap region 
and crossing it, as also seen in their young star distribu- 
tion plot, and presumably extending from the northern 
H I tidal tail to the eastern side of NGC 4039. This may 
indicate that the molecular arm is linked to the base of 
the northern H I tidal tail. A deeper study with fur- 
ther numerical simulations would be needed to match 
the molecular gas morphology found in the observations. 

ALMA allows us to tackle the origin and properties 



of such weaker yet now detectable molecular gas com- 
ponents, which should be widely present in interacting 
systems. Data from other more compact configurations 
of ALMA, including zero and sh ort spacings provid ed by 



the Atacama Compact Array ( jlguchi et al. 
essential to reveal the distribution and kmema 
erties of the gas in such features. 



2009), are 
tic prop- 
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Fig. 1. — HST Optical composite image of Antennae galaxy (Whitmorc ct al. 2010), indicating the positions of the NGC 4038/9 nuclei, 
the overlap region, NGC 4039's spiral arm, and the molecular arm ot tidal origin that we study in this paper. The rectangle indicates the 
mosaic covered by ALMA observations. 
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Fig. 2. — ALMA CO(2-l) moment maps of the southern mosaic. (Upper panel) The moment map, integrated from 1300 to 
1800 km s _1 . We indicate the location of the molecular arm (dashed line), the overlap region and the NGC 4039 nucleus. The boxes 
indicate the regions where we obtained the SF law plot in Figure [3] (Lower panel) Velocity field (moment 1; left), from 1450 to 
1900 km s — 1 ; and the velocity dispersion map (moment 2; right), from to 90 km s — 1 . The beam is indicated as an ellipse in the lower 
left corner of each map (HPBW of l'/68 X 0'.'85 and a PA of 78°). 
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Fig. 3. — (Upper panel) Velocity versus distance along the trace linking the points of higher peak temperature (see lower panel). The 
solid line indicates the centroid velocity of the fit, and the dashed-dotted lines the ± 1/2 dispersion. The grey scale indicate peak intensity. 
The integrated flux along the trace is shown as a dashed profile. The vertical lines indicate profiles where the fit did not converge, given 
initial conditions and constraints to fit. (Lower panel) Integrated intensity map, emphasizing the tidal molecular arm. The contours are 
0.5, 1, 1.5, 2, 3, 5, 10, 22 Jy beam -1 km s" 1 . The plus signs show high local peak temperature as in the upper panel. The (dashed line) 
circles show the region excluded in the SF law plot for the molecular arm region (Section 4 and Figure |4j. The triangle signs show the 
location of super stellar clusters along the molecular arm (see Section 4). We indicate SSC 33, for which a reasonable accurate age could 
be estimated. 
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Fig. 4. — Star Formation Rate (SFR) versus H2 surface densities (Schmidt-Kennicutt law, see Section|3]and|4]l , emphasizing the molecular 
arm (squares), NGC4039 (green circles), and the overlap region (blue circles). Data points for other regions in the southern mosaic are 
shown with plus signs. The rectangles used for each region are indicated in Figure[2) The lines indicate constant SFE, from top to bottom, 
at 10 -7 , 10~ 8 , and 10 -9 yr — 1 . Each data point corresponds to a pixel of 3" in a <?' resolution map. 



